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ABSTRACT 



Very high intensities of galactic cosmic rays measured by Voyager 1 in the 
heliosheath appear to be incompatible with the presence of a modulation "wall" 
near the heliopause produced by a pile up of the heliospheric magnetic field. We 
propose that the modulation wall is a structure permeable to cosmic rays as a 
result of a sectored magnetic field topology compressed by plasma slowdown on 
approach to the heliopause and stretched to high latitudes by latitudinal flows in 
the heliosheath. The tightly folded warped current sheet permits efficient cosmic- 
ray transport in the radial direction via a drift-like mechanism. We show that 
when stochastic variations in the sector widths are taken into account, particle 
transport becomes predominantly diffusive both along and across the magnetic 
sectors. Using a test-particle model for cosmic rays in the heliosheath we inves- 
tigate the dependence of the diffusion coefficients on the properties of the sector 
structure and on particle energy. 

Subject headings: cosmic rays — magnetic fields — turbulence — solar wind 
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Introduction 



The twin Voyager space probes are currently exploring the distant heliosheath region 
and are expected to reach the heliospheric boundary, known as the heliopause, and pass 
into the surrounding interstellar space before the end of the decade. Before crossing the 
heliopause, the Voyagers are expected to enter a region of rapidly changing magnetic field 
which is a remnant of magnetic sectors separated by an oscilla ting heliospheric current 



sheet (HCS). Fr o m glo bal h eliospheric simu 



Borovikov et al. 



(1201 If ), and 



ations reported by 



Czechowski et all (120101 ) 



Florinskil ( 1201 ll ). it emerges that sectors become nar rower as 



the so lar wind plasma slows down as it nears the heliopause. It was suggested by 



Florinski 



( 1201 ll ) that highly energetic charged particles, such as galactic cosmic ray ions, will be 
profoundly affected by the compressed sector structure in the heliosheath. Charged particles 
find it much easier to travel across the sectors, and hence across the magnetic field, when 
the width of a sector is comparable to or smaller than their Larmor radius (rather than 
being much larger than the latter, as is typically the case in the solar wind). The enhanced 
cross-field transport process is a consequence of fast drift-like motion across the stack of 
magnetic sectors. 



Voyager observations (IBurlaga et al.l 120031 ) show no regular magnetic s ector pattern in 



the d istant solar wind, a trend that appears to persist into the heliosheath (IBurlaga fc Ness 
2010). Measured sector durations in the solar wind were, on average, shorter than 13 days 
expected from a simple two-sector pattern implying a significant di s tortio n of the shape of 



the HCS. The widths of magnetic sectors reported by 



Burlaga et al 



(120031 ) varied by almost 



a factor of 10 and had a sufficiently broad distribution. 



Burlaga et al. 



( 120031 ) only reported sectors observed during three years near the 
maximum of solar activity. Close to solar minima magnetic sectors occupy only a narrow 
range of heliolatitudes. They were not observed by the Voyagers during the pre-termination 
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shock part of the mission when the spacecraft was above (VI) or below (V2) the current 



sheet. The situation changes dr amatically in the 



away from the radial direction ([Richardson et al. 



and longitude (ICzechowski et al. 



2010; 



leliosheath. Here the flow is deflected 
2009|) stret ching the HCS in latitude 



Borovikov et al. 



20111 ) . In the hemisphere facing 



the interstellar flow the sector structure is expected to exist at all latitudes deep into the 
heliosheath. It is then of great importance to understand the effect of the tightly folded 
current sheet on energetic particles, such as galactic cosmic rays. 



It was proposed in iFlorinskil ( 1201 ll ) that variability in sector width results in an 
enhanced cross-field diffusion of galactic cosmic rays with a Larmor radius greater than 
about 0.1 AU. The process could well be responsible for very high intensities of galactic 
helium ions measured by Voyager 1, inten sities that are presently ( 2011-2012 ) near or in 



excess of the predicted interstellar values ( j Webber fc Higbie 



20091 ). 



Florinskil ( 1201 ll ) used 



a simple, one-dimensional sector model that did not take into account variations in the 
shape of current sheets. The present paper reports a much more complete, two-dimensional 
model of particle transport near a rapidly oscillating HCS. Here we assume that the sector 
width in the distant heliosheath varies about a well defined mean, and that the variation 

acement in the 



Florinski 



( 1201 lh 



occurs as a result of differential (variable in latitude) current sheet disp 
radial direction. This work therefore builds upon the ideas proposed in 
to develop a quantitative model of cosmic-ray transport in a sectored magnetic field of the 
distant heliosheath. 

It is noteworthy that a partially related concept was propos ed much earlier as a 



mechanism of interstellar dust grain transport in the solar wind (jMorfill fc Griin 



Czechowski fc Mann 



1979 



20031 ). High rigidity and small velocity of dust grains makes their 



transport qualitatively different from that of cosmic rays. Dust particles larger than about 
0.01 fim have gyroradii comparable to the size of the heliosphere and their trajectories 
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are only mildly perturbed by the sector structure in the solar wind. These slowly moving 
particles acquire a velocity component perpendicular to the interplanetary magnetic field 
and parallel to the sectors as a result of the motional electric field in the frame of the 
particle. Random variations in sector duration then lead to a diffusion in pitch angle. For 
galactic cosmic rays with their high velocities this effect is negligible, and the principal 
mode of transport is drift associated with the current sheet crossings. 

The paper is structured as follows. In Section 2 we discuss the fundamental physics of 
particle drift in a sectored magnetic field. In section 3 we investigate cosmic ray transport 
through a stack of magnetic sectors with a test-particle orbit integration method commonly 
used to simulate scattering and cross- field transport in tur bulent magnetic fields (e.g., 



Giacalone fc Jokipii 



1999 



Mace et al. 



2000 



Qin et al. 



20021 ) . We demonstrate that particle 



transport in a magnetic field composed of variable-width sectors of alternating polarity is 
diffusive both parallel to and across the current sheets. The results from this study are 
reported in Section 4. These results reveal the dependence of the diffusion coefficients on 
the distribution of sector widths and on the cosmic-ray energy. Finally, Section 5 explores 
some consequences of the new diffusion mechanism for cosmic-ray propagation through the 
heliosheath. 



2. Neutral sheet diffusion as a drift process 

Consider an energetic charged particle, such as a galactic cosmic ray ion placed in a 
magnetic field that has the same magnitude everywhere, but can change sign. A change 
of sign occurs across a neutral sheet with negligible thickness (i.e., thickness much smaller 
than the Larmor radius of the energetic particle). It is well known that in the absence of a 
large-scale plasma flow, the particle can either gyrate in place (ignoring the motion along 
the magnetic field, which is of no interest here), or undergo a meandering motion (drift) 
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along a single current sheet. When the sectors become squeezed together, as one expects 
will happen on approach to the heliopause, the distance between adjacent neutral sheet 
folds would become comparable or smaller than the cosmic ray cyclotron radius r g , which 
is typically between 0.05 and 0.3 AU in the heliosheath. 

The slowly moving heliosheath plasma is th ought to be nearly incomp ressible, which 



2011 



). The plasma is 



appears to be consistent with Voyager 2 results ([Richardson &: Wangi 
therefore not compressed, but follows different paths toward the tail (spreads out in latitude 
and longitude). Evidently, for primarily north-south flows the magnitude of the mean field 
need not increase because the decrease in the radial distance between field lines (as the 
sectors become thinner) is balanced by spreading of the field lines in latitude. For azimuthal 
flows, the magnitude of B should increase because in that case plasma follows magnetic 
field lines. Altogether, we expect the magnetic field to increase by a modest amount, so 
that the ratio of r g to the mean distance between sectors d increases with heliocentric 
distance. At some distance, therefore, the critical condition 2r g = d will be satisfied. 

To support our estimate we calculated the magnetic field and the mean sector width 

in the heliosheath as a function of radial distance from the Sun using a three-dimensional 

I II 1 

model of the heliosphere that includes neutral hydrogen atoms (jFlorinskil 1201 ll ). From 
this data we calculated the Larmor radius r g for 20, 200, and 2000 MeV/n He 2+ ions as 
a function of heliocentric distance. Figure [1] shows the radial dependence of d and 2r g in 
the upwind (nose) direction. Because of its blunt shape, the termination shock stand-off 
distance, marked by the dotted line, is smaller (~ 74 AU) than in the directions of Voyager 
1 and 2. The model is steady-state and does not include dynamic solar wind structures, but 
is sufficient to illustrate our point. Readers interested in the mo del boundary co nditions for 



the solar wind and the local interstellar medium are referred to 



Florinskil ( 1201 lh 



From the figure, the gyro-radius does decrease with distance, but not as rapidly as 
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Fig. 1. — Larmor radius of galactic helium at three different energies per nucleon (solid 
lines) and the mean thickness of a magnetic sector in the heliosheath in the upwind (nose) 
direction. Current sheet diffusion is expected to be important between the intersection point 
of the 2r g and d curves and the heliopause. The positions of the termination shock and the 
heliopause are marked with dotted lines. 
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the sector width. In the model the magnetic field magnitude increases from 0.15 nT to 0.5 
nT between the termination shock and the heliopause, whereas the sector width decreases 
by nearly two orders of magnitude. The point of intersection between the solid and the 
dashed curves in Figure [I] {2r g = d) gives the approximate distance beyond which the 
effects of neutral shift diffusion become noticeable. For 200 MeV/n ions this occurs about 
two-thirds of the way through, but only close to the heliopause for 20 MeV/n and lower 
energy particles. 



As discussed in iFlorinski 



(l201lf ). for r g < d/2 there exist gyrophase angles for which 
meandering or gyrating particles are unable to complete their cycle before encountering 
a neighboring current sheet. At that point the particle's sense of rotation will change as 
its guiding center "jumps" across the sheet and into an adjacent sector. Figure |2] shows a 
typical trajectory of an energetic particle in an ensemble of sectors of uniform width. The 
particle drifts in both y and x directions. However, unlike the former drift, it is easy to 
see by drawing a similar diagram for a different starting gyrophase, that drift across the 
magnetic boundaries can proceed in either direction (to the left or to the right on the figure, 
depending on the gyrophase). It follows that the net drift velocity is zero for a gyrotropic 
ensemble of particles. However, the process itself may lead to a separation of an initially 
symmetric distribution into left- and right-going components. 

Despite the limitation, the drift process discussed above could enhance the diffusive 
transport of energetic particles if the sector structure is not uniform. Observations clearly 
show that the sector widths (or duration, in terms of spacecraft measurements) vary 



by a factor of ~ 10 in the outer heliosphere (IBurlaga et al. 



20031 1 . This variability is 



presumably due to both an irregular magnetic field structure on the Sun and the dynamic 
interaction between structures in the sola r wind further distorting the shape of the HCS (e.g. 



Riley et al 



2002 



Intriligator et al. 



20051 ). One expects the sector thickness to vary with 
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Fig. 2. — Drift motion of a positively charged ion through a stack of magnetic sectors of 
alternating polarity. The solid line is the particle's trajectory and the dotted line shows the 
"trajectory" of its guiding center. The magnetic field is in and out of the plane of the figure. 
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both latitude and longitude; we include only latitudinal variations in our two-dimensional 
model (see below). The advantage of the present model is simplicity stemming from the 
assumption that the magnetic field varies in direction only, but not in magnitude (this 
condition cannot be maintained in the fully three-dimensional case). This allows us to focus 
on the fundamental physics of the diffusion process due to current sheet fluctuations, a 
feature that would otherwise be hidden underneath the complexity of a fully 3D simulation. 

In the last part of this section we present a simple two-dimensional model, based on 
the Boltzmann equation, that illustrates the qualitative features expected from closely 
packed alternating sectors. In this model, the sector width is fixed and the randomness is 
provided by scattering centers embedded in the plasma. For the sake of brevity, we consider 
only the motion in the xy plane (with the background magnetic field directed along z), 
and disregard scattering in pitch angle, which would lead to quantitative modifications but 
would not change the robust qualitative features. We assume the pitch angle to be 7r/2 and 
consider scattering in gyrophase angle ip only. 

As in Figure [2] we take the distribution function of cosmic rays to be independent 
of y. Assuming that diffusion is a result of independent small-angle scattering in if, the 
Boltzmann equation for the phase space density f(x,ip,t) can be written as 

dt V C ° S ^ dx dip dip \t dip J ' 

where Q = v/r g is the gyro- frequency and r is the mean scattering time. The velocity v 
is the same for all particles. Physically, r _1 is related to the quasi-linear Fokker-Planck 
coefficient D w that descr ibes diffusion in gy rophase through resonant interactions with 



turbulent magnetic fields (lAchatz et al. 



19911 ). The latter is a function of pitch-angle cosine 



\i and velocity v; since both are constant in our model, we accordingly take r = const. 

For a unipolar field pointing in the z direction (i.e. constant Q), Eq. (1) entertains the 
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following steady-state solution 



25*° 2S° 
f(x, ip) = F + G x x H cos ip H sin <p, 



(2) 



implying that a steady gradient (G x ) produces streaming S^. y = (2n)~ 1 j fv x<y d(p = 
—H( x ,y)xG x in the x and y directions, respectively, consistently with the nominal diffusion 
coefficients of k xx = v 2 t/[2(1 + Q 2 t 2 )} and k 



-VtTK n 



Next, we place alternating sectors of widths d + and oL periodically in the yz plane. 
Though the solution shown in Eq. (2) remains formally a solution, one should bear in mind 
that S y now becomes discontinuous at the current sheets, while the p hysical solutio n must 



19651 ) cannot 



be continuous. We note that an approach based on Parker's equation (IParkerl 
be directly applied near current sheets because Parker's equation does not include shear. 
Its extensio n to include second-order te rms shows that shear produces significant higher 



harmonics (IKota 



1975 



Earl et al. 



19881 ) . which cannot be neglected in our case. 



To represent a large scale gradient, G x , we assume a solution of the Boltzmann equation 
(1) in the form: 

f(x,(p) = G x x + F(x,(p), (3) 

where F(x, if) is periodic in x and is required to be continuous at the currents sheets. 
Integrating Eq. (1) immediately yields a physically plausible result that streaming in the x 
direction, S x , must be constant. The ratio of the computed S x /G x represents the effective 
diffusion coefficient appropriate to this situation. We solve Eq. (1) numerically on a 
200 x 120 grid in x and if. Numerical results presented here were all obtained with A = 5r g 
(i.e., Qt = 5). Furthermore, we consider a symmetric arrangement with d+ = cL = D/2, as 
well as an asymmetric one with a dominant sector by choosing d + = 3D/4 and gL = D/A. 

Figure [3] shows how the streaming S x changes if the sectors become narrower, i.e., 
r g /D increases. Because S x ~ k xx (assuming a constant spatial gradient), the figure also 



- 12 - 




Fig. 3. — The increase of the streaming S x in the x direction relative to the nominal streaming 
obtained for a uniform field. The solid curve refers to equidistant alternating sectors, the 
dashed curve refers to an asymmetric case with a dominant sector (see text). 
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Fig. 4. — Spatial variation of the diffusive streaming S y (solid: r g /D = 0.02, dotted: r g /D = 
0.2, dashed: r g /D = 2). While S y (x) has a sharp transition at the HCS for r g /D 1, it 
decreases in magnitude as r g /D increases. The left panel shows a symmetric case with two 
equally wide sectors, the right panel shows the asymmetric arrangement with a dominant 
sector (see text). Note that S y remains large in the latter case. 
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shows the ratio between the cross-field diffusion coefficients with and without the sectors. 
The two curves refer to the symmetric and asymmetric arrangements, respectively. Clearly, 
particle transport is significantly faster as the thickness of the sectors becomes comparable 
or smaller than the particle's gyro-radius, r g . The increase is larger in the symmetric case, 
when S x approaches the limiting case obtained without any background field (Q = 0), 
which is about a 25- fold increase for our parameters (Qt = 5). For the asymmetric case 
the increase is slower. Intuitively one can see that particles will sense the average Q, which 
implies a twofold decrease in (Q), and a corresponding fourfold increase in S x relative to 
the unipolar case. 

Figure H] illustrates the spatial variation of the S y streaming component for the values 
of T g jD = 0.02, r g /D = 0.2, and r g /D = 2 within a single pair of sectors (the complete 
solution is of course periodic in x). For small r g /D the solution is close to the formal 
discontinuous solution, with a sharp transition occurring near the current sheet. We note 
that the directional distribution has significant higher order anisotropies in this region. S y 
decreases as r g /D increases. The right panel of Figure H] shows the asymmetric case, with 
one sector dominating. This case yields a remarkably large streaming in the y-direction 
in both sectors for r g /D = 0.2. This seemingly counter-intuitive result can be readily 
understood recalling that multiplying Eq. (1) by sin tp and averaging over x gives 

(S y ) = (nr)S x , (4) 

which means that large streaming in the x direction implies large streaming in the y 
direction as well, if one sector dominates. 
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3. Test particle simulations of neutral sheet diffusion 



Building upon the conceptual random sector thickness model introduced in iFlorinski 



( 120111 ). we now proceed to investigate the transport of energetic particles across a more 
realistic representation of magnetic sectors in the heliosheath using a test-particle computer 
code for orbit integra tion. Measured plasma velocities are small in the distant heliosheath 



( jKrimigis et al. 



201 ll ). and we ignore the background flow altogether. Inside the sectors 
we assume that the magnetic field B is constant in magnitude and points either in the +z 
or the — z directions, corresponding to the primarily azimuthal heliosheath magnetic field. 
The motion of a charged particle is unrestricted along the field, and we only consider their 
trajectories in the xy plane (in other words, z is an ignorable coordinate in the present 
model). 

In the computer model used here particles are injected into a two-dimensional 
rectangular simulations box of size L x x L y containing a few hundred alternating polarity 
sectors separated by current sheets (which represent different segments of the same HCS). 
Whereas a theoretical number of sectors in one 11-year solar cycle is about 150, we use a 
larger number in the model to obtain a correct asymptotic behavior for large times. The 
current sheets are, on average, parallel to the yz plane, and their mean spacing along the 
x (radial) direction is d. We introduce perturbations to the form (x-displacement) of each 
current sheet according to 



N k 



x s (y) = x s o + ^2 A n cos(k n y + (3 sn ), 



(5) 



n=l 



where x s q = sd (s is the integer label enumerating the sheets), k n is the wavenumber of 
mode n, (3 sn is the random pha se, and A„ is the amplitud e, related to the power spectrum 



of the fluctuations S(k) as (e.g. 



Giacalone &: Jokipii 



1999|) 



A 



S{k n )Ak n a 2 x 



(6) 
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The model one- dimensional spectrum of current sheet displacements consists of a flat 
energy range between fc min and the "bend-over" wave number, followed by an inertial 
range with 5* ~ (k/kb)~ 5 ^ 3 extending to fc max - We use N k = 150 logarithmically equidistant 
modes separated by Ak n . For the heliosheath we assume that A; max = 2vr/0.05 AIT 1 and 
kmm = 27r/79 AlH 1 , a difference of 3.2 orders of magnitude, with kb at the midpoint of the 
logarithmic interval. This is our "standard" sector boundary displacement structure. From 
the central limit theorem one may infer that the sector width distribution thus obtained 
should be Gaussian, with mean d and dispersion a x . A direct numerical calculation of the 
distribution function confirmed that this is indeed the case. 

In our simulations length is measured in units of the largest (for pitch angle vr/2) 
particle gyro-radius r g = pc/(ZeB), where p is momentum and Z is its charge number. 
The field is computed on a fine grid with equal mesh spacing in both directions; a typical 
simulation may contain up to 30 billion grid cells. The L x : L y aspect ratio varies depending 
on the ratio of calculated diffusion coefficients in the two directions, k xx and ^>yy ^SGG 
below). Typically we set L x /L y {k xx / K yy ) 1 / 2 . The size of the box is chosen such that no 
particles escape through any side of the box by the end of the run time. Note that particles 
traveling the farthest from their starting point contribute the most to diffusive transport 
statistics at later times. If such particles were allowed to leave the system, the resulting 
diffusion coefficients would be underestimated. Figure [5] shows a fragment of the sector 
boundary structure calculated with d = 1.25r 9 and a x = 0.3r 9 . The wiggles in the shape of 
the sheet are assumed to be a result of turbulent flows in the solar wind and heliosheath 
that displace the HCS radially. 

In our model the sector structure is controlled by three parameters: d, a x , and kb 
(assuming the ratio k mauX /k min is a constant). The resulting magnetic field topology may be 
viewed as a special case of strong "turbulence" with zero mean field. This type of turbulence 
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Fig. 5. — Magnetic sector boundaries (a small fragment of the entire simulation) in the xy 
plane with d = 1.25r g , a x = 0.3r g , and k b = \.h7r~ l . The magnetic field is perpendicular to 
the page. 
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does not belong to the conventional slab/2D classification system ( IMatthaeus et al.lll990 



Bieberetal. 



19961 ). In our case the fluctuating field <5B || z is orthogonal to the plane 
spanned by the wavevectors (the xy plane). The turbulence is clearly non-axisymmetric; 
the perpendicular (x) power spectrum is dominated by a signal at ir/d (twice the mean 
current sheet spacing), and the parallel (y) spectrum has no clearly identifiable dominant 
wavenumber. 

Charged particle orbits are computed by integrating the Newton-Lorentz equation of 
motion 

f = vxn, (7) 

where v is the particle's velocity and fl = ZeQ / (7777,0) is its cyclotron frequency, m being 
the mass and 7 the relativistic factor. Eq. (7) is solved numerically using the standard 
fourth-order Runge-Kutta (RK4) method. To achieve good accuracy we used a small time 
step, typically At = 27r(1000fi)~ 1 . The spatial resolution in both (x, y) directions was 
Ax = Ay = 0.025r 5 (dependent on particle's energy). 

Figure illustrates typical particle behavior in a sectored magnetic field with 
stochastically displaced boundaries. The particle was injected with a fixed pitch angle 
9 = 90° near the center of the box. Notice the long excursions in the y direction, 
corresponding to particle drifting along a single current sheet for long periods of time. The 
particle may subsequently transfer to an adjacent sheet when the two become sufficiently 
close. By examining Figure [6] one can envision a certain parallel between particle motion 
along the current sheets vs. motion across the sectors and motion along and across the 
mean magnetic field in a conventional slab/2D composite turbulence. Indeed, in the latter 
case a particle stays on the same field line for a long time, but can at times switch to a 
different field line as a result of pitch-angle scattering events randomizing its gyrophase. A 
constriction in the sector structure acts similar to a scattering center allowing a particle to 
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transfer from one current sheet to another. 

From our simulations it emerges that particle transport in both directions becomes 
diffusive after a certain period of time (dependent on particle's energy, mean sector width, 
and sector displacement properties). We use running diffusion coefficients calculated as 

^ (t) = <WQ^, ^ t) = Mg£), (8) 

where (xo,yo) are the initial coordinates of the particles. Here angular brackets mean 
ensemble averaging over both particles and field realizations. We typically use 25 sector 
structure realizations and 2000 particles per realization injected in an isotropic shell 
distribution with a specified kinetic energy T. Figure [7] shows k xx and K yy as functions of 
time for a demonstration run for isotropically injected 200 MeV/n He 2+ in a field with a 
strength B = 0.4 nT, and the sheet structure as in Fig. |5j The size of the simulation box 
was 107 AU by 573 AU in the x and y directions, respectively. Such a large size was chosen 
to ensure that we measure correct asymptotic values of the running diffusion coefficients 
(see discussion above). All 50,000 particles remained inside the box after 2000 gyro-periods. 

After a brief super diffusive phase, lasting a few tens of gyro-periods, when particles 
are streaming primarily along individual neutral sheets (marked by increasing n(t)), normal 
diffusion is established in both directions. Asymptotic values of k xx and obtained 
from this late phase of the particles' trajectories. For the present choice of parameters, 
the cross-sector diffusion coefficient is a factor of ~ 4 smaller than the "parallel" diffusion 
coefficient. The ratio K xx /n yy depends sensitively on the width of the sectors. The next 
section demonstrates the dependence of diffusion coefficients on particle's energy and on 
the properties of current-sheet "turbulence" . 
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Fig. 7. — Running diffusion coefficients, in units of r^Q, as a function of time. Isotropically 
injected 200 MeV/n helium ions in a 4 /iG magnetic field, mean sector width d = 0.062 AU. 
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4. Dependence on sector properties and particle energy 

We first discuss the dependence of diffusion coefficients on energy. In the following 
numerical experiments the current sheets had the same properties using the "standard" 
choice of parameters B, d, a x , k min , k b , and fc max as described in the previous section. The 
energy of He 2+ ions was varied between 10 MeV/n and 1 GeV/n. The size of the simulation 
box was adjusted accordingly to prevent particles from leaving. 

The results of this experiment are shown in Figure El At the lowest energy (10 MeV/n) 
the Larmor radius is about 4 times smaller than the mean distance between the sheets 
and perpendicular diffusion is strongly suppressed. Parallel transport still appears to be 
super-diffusive by the end of the simulation, which may explain the upturn in n yy at low 
energies. By about 100 MeV/n the neutral sheet separation is about equal the Larmor 
radius leading to a dramatic increase in k xx . Here transport is strictly diffusive and both 
coefficients increase with ion energy. The parallel diffusion coefficient is well approximated 
by a power-law with an index of 3/2, while the perpendicular coefficient increases somewhat 
faster as ~ T 2 . In the non-relativistic regime, this behavior is comparable with a q uasi-line ar 
theor y result n yy ~ j^W 2 f or a turbulent power-law spectral slope a = 1 (e. 



Jokipii 



197ll ). It is not clear whether there is a connection between the two or it is a simple 
coincidence. Note that in our model the spectrum of fluctuations in the current sheets 
varies between —5/3 and zero. 

Figure [9] shows the dependence of k xx / k w on ion energy. The ratio is a rapidly 
increasing function of energy (~ T 3 at low energy) and approaches a unity (isotropic 
diffusion) at high energy. The transition from slow to fast perpendicular diffusion occurs at 
r g ~ d. 

In the second set of numerical experiments we simulated particles at two fixed 
energies, 50 MeV/n and 250 MeV/n, while varying the distance d between the sheets and 
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Fig. 8. — Parallel (/%,) and perpendicular (k xx ) diffusion coefficients for galactic He as a 
function of energy per nucleon. The dashed line has a power law slope of 3/2. 
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T, MeV/n 

Fig. 9. — Ratio of perpendicular (k xx ) to parallel (K yy ) diffusion coefficients for galactic He 
as a function of energy per nucleon. The dashed line has a slope of 3. 
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the displacement magnitude a x while keeping the ratio o x jd a constant (= 0.24). This 
approximates a situation where the stack of sectors is compressed, but the wiggles become 
smaller as well (the complete structure is flattened on approach to the stagnation point). 
The magnitude of B was kept a constant implying that the excess field was removed from 
the region by a latitudinal flow. The current sheet turbulence spectrum was again of the 
"standard" type (described following Eq. (6)). 

We show the ratio K xx /n yy in Figure [TU1 For 50 MeV/n ions the ratio drops rapidly as 
the distance d becomes larger than the gyro-radius. Conversely, at 250 MeV/n the decrease 
is much more gradual. The main reason for this difference is the fact that the power 
spectrum of neutral sheet oscillations is the same in both simulations (with 50 MeV/n and 
250 MeV/n). The two would be alike only if the neutral sheet displacement wavelengths 
were also varied ~ d. 

In the final experiment we vary the turbulence spectrum for a fixed energy (100 
MeV/n) and sector width (d = 1.25r 9 ). Instead of a broadband spectrum employed in the 
previous simulations, here we use a narrow-band spectrum with a constant power spectral 
density between k min and fc max . The spectrum extends a single decade in wavenumber and 
is represented with iV& = 50 modes. Running diffusion coefficients at the end of each run 
are recorded as the asymptotic values. 

The results are plotted in Figure [TTJ We see that increasing the wavelength leads to an 
increase in K yy and a decrease in k xx . The latter is easily understood from the fact that by 
reducing h m ™ the mean distance between neutral sheet constrictions, acting as scattering 
centers, increases leading to a reduction in scattering rates. There is a rapid drop in K yy 
near what could be described as "resonant" wavenumber where kr g is of order 1. It is not 
clear what causes the upturn at large /c m j n . It is conceivable that the model is unable to 
resolve such fine structures with good accuracy. 
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Fig. 10. — Ratio of perpendicular (k xx ) to parallel (/%,) diffusion coefficients at 50 MeV/n 
(squares) and 250 MeV/n (diamonds) as a function of mean sector thickness. The ratio of 
the magnitude of sector boundary fluctuations and the sector thickness was kept a constant 
(0.24). Vertical dashed lines are the gyro-radii of both ions in a 0.4 nT magnetic field. 
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Fig. 11. — Parallel (/%,, diamonds) and perpendicular (k xx , squares) diffusion coefficients of 
100 MeV/n He ions as a function of k min . The wavenumber and the diffusion coefficients are 
normalized to r g and r^Q, respectively. The mean sector width d = 1.25r g . 
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5. Discussion 

Analytic and numerical results presented here clearly demonstrate that a warped HCS 
can facilitate charged particle transport perpendicular to the mean magnetic field, if the 
distance between the folds d is comparable to (or smaller than) the particle's Larmor radius, 
the condition that is expected to be satisfied within about 10 AU before the heliopause. 
The analytic treatment of the problem was performed using the Boltzmann equation for 
a non-gyrotropic velocity distribution. Unlike the numerical model, the sector width was 
constant. Instead, stochasticity was provided by small-angle gyrophase scattering events. 
Our results show a steady increase in the cosmic-ray differential current perpendicular 
to the sectors (and therefore, in the cross-field diffusion coefficient) with the ratio r g /d. 
Gyration can produce net streaming along the sector planes if the widths of the positive and 
the negative sectors are different. The magnitude of the streaming flux tends to increase 
with r g /d. 

In our test-particle simulations we found that particles undergo diffusive random 
walk in the plane perpendicular to B, i.e., in the radial and latitudinal directions. The 
efficiency of this process depends critically on the size of neutral sheet "wiggles" and the 
mean thickness of the alternating polarity sectors. The ratio of perpendicular and parallel 
diffusion coefficients is zero in the limit of very long-wavelenfth fluctuations, and when 
the ratio r g /d is small. In the opposite limit (short-wavelength fluctuations, r g /d > 1) 
the ratio approaches unity. In the latter regime, the value of K rr could be as large as 
6 x 10 21 — 6 x 10 23 cm 2 /s for He 2+ ions with energies between 100 MeV and 1 GeV per 
nucleon, assuming a magnetic field strength of 0.4 nT (see Figure [H])- These values are 
larger than the perpendicular diffusion coefficients calcula ted using the standard theor y 
for solar-wind fluctuations. For example, in the model of 



Florinski fc Pogorelov 



([20091 ) 



k± is between 4 x 10 21 and 4 x 10 22 cm 2 /s for 100 MeV and 1 GeV protons, respectively. 
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For helium with the same energy per nucleon the diffusion coefficients will be larger by 
up to a factor of 4 (for resonant interactions in the energy range of the turbulence with 
K\\ ~ r 2 g and k± ~ Ky). The current-sheet radial diffusion coefficient is about a factor of 5 
larger. Remember, however, that conventional scattering by turbulence was not included 
in our model. We expect our results to be more reliable in a situation where the level of 
background turbulent activity at resonant wavenumbers is relatively low. 

The interplay between the sheet-parallel and perpendicular transport in our model 
resembles in several ways conventional diffusion along and across the magnetic field. 
Constrictions between the adjacent sheets appear to play the role of scattering centers. 
Nevertheless, one should exercise caution before drawing parallels between the two. Parallel 
diffusion in QLT is a resonant process, which is clearly not the case for current sheet 
diffusion. 

The current sheet transport mechanism is most efficient for galactic ions with energies 
of the order 100 MeV per nucleon and above. Voyager 1 observations of galactic helium 
made near the end of 2011 at a heliocentric distance of ~ 119 AU (McD onald et al., AGU 



abstra ct) are close or in excess of the predicted interstellar intensities from 



Webber fc Higbie 



( 120091 ). If Voyager 1 is 10 AU in front of the heliopause, the current-sheet effect could be 



responsible f or the lack of modu 



ation. Another possible cause is the absence of a radial 



20111 ) which would otherwise expel energetic particles from the 



plasma flow ( IKrimigis et al. 
heliosheath. 

The diffusion mechanism discussed here may facilitate the transport of interstellar dust 
grains with size < 0.01 Atm that have a gy roradius similar to that of a few hundred MeV 



cosmic-ray ion (ICzechowski fc Mann 



20031 ). Conversely, it is not efficient for low-rigidity 
particles such as electrons, because their gyro-radius is much smaller than the distance 
between the adjacent sheets. The recent rise of the low-energy electrons measured by the 
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Voyagers (jCaballero-Lopez et al.ll2010l ; IWebber et al.ll2012l ). is probably not directly related 



to the current sheet diffusion effect discussed here. At such small rigidities the effect of 
magnetic field meandering is likely to take precedence over current sheet transport as 
described by our current model. The model is somewhat restricted by the assumption that 
the direction along B is an ignorable coordinate; this could lead to an under-estimate of the 
diffusion rate in the radial direction. 

We realize that a more accurate diffusion model should include "conventional" 
turbulence (variations in the magnitude and direction of B) in addition to the current-sheet 
oscillations. Such a project would require massively larger computer simulations taking 
thousands of CPU-hours for a single run, but may be ultimately needed solve the problem of 
cosmic-ray transport in the distant heliosheath. Despite technical difficulties, in the future 
it might be possible to combine the two kinds of diffusion into a single three-dimensional 
model of energetic charged particle transport. 
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